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ABSTRACT 



A optical analytical method and apparatus for measuring 
concentrations of materials in liquids and gases is described. 
The method is based upon measurements of light reflection 
from diffraction gratings in contact with the sample. Fea- 
tures in the reflection spectrum are directly related to the 
complex dielectric constant of the sample solution. By 
examining the reflection spectrum, the concentration of 
materials in the liquid or gas can be deterrnined. 



11 Claims, 7 Drawing Sheets 
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ANALYTICAL SENSOR USING GRATING 
LIGHT REFLECTION SPECTROSCOPY 

FIELD OF THE INVENTION 

The present invention relates to analytical sensors and, 5 
more particularly, a sensor that uses grating light reflection 
spectroscopy to determine the concentration of a material in 
liquids or gases. 

BACKGROUND OF THE INVENTION 10 

In recently published articles and patents there were 
described a number of methods for measuring bulk fluid 
compositions in statis and flow by visible light reflection 
characteristics from the optically inhomogeneous wall of a 15 
fluid containing vessel. Surface Plasmon Spectroscopy 
(SPS) on gratings and randomly rough surfaces has proven 
to be especially sensitive to fluid composition near the wall. 
See M. J. Jory, P. S. Vukusic, J. R. Sambles, Sensors and 
Actuators B: Chemical 17, 203-209 (1994). It has been 2 o 
shown that SPS combined with electrochemical modulation 
techniques allow for the determination of relative concen- 
trations of surface active species as low as 10~ 15 to 10~ 14 M. 
See A. M. Brodsky, M. Urbach,&>v. Phys. Uspechi, 132, 413 
(1982). 25 

SPS techniques in general have a number of inherent 
drawbacks when used as analytical sensing techniques. For 
instance, the metal dielectric interface that is necessary for 
surface plasmon generation must be chosen carefully. As 
only a limited number of metals (gold, silver) are available 30 
for the SPS, the technique is limited to environments that 
will not corrode or destroy the surface plasmon generation 
medium. In contrast, the technique described herein, termed 
Grating Light Reflection Spectroscopy (GLRS), relies on an 
optical phenomenon that requires only a sufficient dielectric 35 
difference between the substrate and "sensing" layers. Thus, 
any combination of substrate and grating materials with 
suflicient dielectric difference will suffice, allowing for 
sensing in environments that would preclude the use of SPS. 

Additionally, the surface sensitivity associated with SPS 40 
techniques make them unstable for real world analytical 
applications as the bulk determinations of solutions would 
be unattainable in the presence of surface active analytes. 
The techniques described herein proves to be relatively 
insensitive to surface active species. Although, here, we 45 
have not shown this effect experimentally, the mathematics 
described below indicate a relative surface insensitivity. 

SUMMARY OF THE INVENTION 

A optical analytical method and apparatus for measuring 50 
concentrations of materials in liquids and gases is described. 
The method is based upon measurements of light reflection 
from diffraction gratings in contact with the sample. Fea- 
tures in the reflection spectrum are directly related to the 
complex dielectric constant of the sample solution. By 55 
examining the reflection spectrum, the concentration of 
materials in the liquid or gas can be determined. 

In particular, the intensity and phase modulations of the 
reflection spectrum are examined to determine the bulk 60 
dielectric properties of the solution or gas in contact with the 
diffraction grating. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing aspects and many of the attendant advan- 65 
tages of this invention will become more readily appreciated 
as the same becomes better understood by reference to the 



,560 

2 

following detailed description, when taken in conjunction 
with the accompanying drawings, wherein: 

FIG. 1 is a schematic diagram of grating light reflection 
spectroscopy (GLRS) on a planar wall; 

FIG. 1A is an illustration of an apparatus formed in 
accordance with the present invention; 

FIG. 2 is a graph of response reflection coefficient for TM 
polarization during grating reflection spectroscopy of EtOH 
solutions in water; 

FIG. 3 is a graph of response reflection coefficient for TE 
polarization during grating reflection spectroscopy of EtOH 
solutions in water; 

FIG. 4 is a graph of PLS prediction of refractive indices 
of EtOH solutions in water for TM polarization in a four 
factor model; 

FIG. 5 is a graph of PLS prediction of refractive indices 
of EtOH solutions in water for TE polarization in a five 
factor model; and 

FIG. 6 is a graphical comparison of the theoretical posi- 
tion of a singularity with experimental spectrum for GLRS, 
TM polarization, and 30% EtOH solution. 



DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENT 

The detailed description is organized as follows. In Sec- 
tion I the theoretical underpinnings of GLRS is described. In 
Section II, the experimental device and procedures is 
described. In Section HI, the experimental results and analy- 
sis for ethanol- water mixtures is presented. In Section IV, 
possible applications of the present invention is described. 

I. Theoretical 

The present invention relies on the propagation modula- 
tion of a light beam that is incident on a diffraction grating 
101 fabricated of metal lines on a dielectric 103, in this case 
chrome and glass. Certain features in the diffracted reflection 
spectra, namely intensity and phase modulations, are 
directly related to bulk dielectric properties of the solution in 
105 contact with the grating (see FIG. 1 ). The external 
media 107 and the dielectric 103 are characterized by 
dielectric functions e (0) (co)and e (1) ((u), respectively, where © 
in the arguments designates the frequency dependence of the 
dielectric function. As the external medium for these experi- 
ments is air, e (0) (co) is taken to be unity for this theoretical 
development. The solution of 105 is characterized by a 
dielectric function e (2) (co) and supports different transmitted 
beams. The emergence of reflected and transmitted beams 
with m*0, where m is the diffraction order 0,±1,±2 . . . , is 
the result of diffraction from the grating structure on the 
bottom of the substrate and in contact with the solution. For 
the sake of notational consistency, the superscripts on the 
wavevector or dielectric function symbols in the following 
equations refers to the medium, as designated above by (0), 
(1), (2), with (0) being the external media 107, (1) being the 
substrate 103, and (2) being the solution 105. Subscripts in 
the equations refer to the cartesian coordinate component, 
with 1 being the normal to the incident plane (see FIG. 1 ). 

The wave vectors k (0 (i=0,l,2) in different media have 

normal k and parallel k 2 (0 components to the surfaces 
which are connected by the following condition: 

c 2 
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where c is the speed of light. If the grating is strictly periodic 

then the parallel components k 2 (i) (m) of the wave vectors 
corresponding to different diffraction orders, m, due to 
momentum conservation laws, are subject to the following 
constraint: 



lU, (m) = -T- sine + m-^- ;m = 0,±l,±2 



(2) 



where a is the grating period, 0 is the incidence angle, and 
m in the argument of the wavevector denotes the difference 10 
wavevector associated with each diffracted order. As it 
follows from equations (1) and (2) the normal components 

k J 2 \m) of the wave vector in medium 2 are equal to: 



*P(m) = 



(3) 



If the value of Ree^o), (Re and Im denoting the real and 
imaginary parts of the complex dielectric function respec- 
tively), is changing composition of the solution of interest 
(medium 2) or due to changing frequency co there exists for 
some value m cn a critical diffraction order, a threshold 
where the function 



(4) 



passes through zero and becomes negative with further 
decline of Re(e (2) (co)). Here, X is the free space wavelength 
of the incident radiation. For the case of 30 
ImCe^^ycReCe* 2 ^©)) we will observe specific singulari- 
ties in the behavior of reflected light The important point is 
that the positions of such singularities, 5m cr =0, and the 
functional dependence on 8m cr in the internal where Sm cr is 
small, I8mj«l, are dependent only on the bulk properties 35 
of the solution at distances from the solution-grating inter- 
face on the order of a few wavelengths. 

The physical origin of such singularities can be inter- 
preted as follows. The intensity I( r ,m) of the light corre- 
sponding to the mth diffraction order in medium (2) is 40 
proportional to: 



/(r,mHE r2J f7,m) E w fT.m)) for r t >0 



(5) 



45 



r is the direction vector for the propagation order as defined 
by the coordinate system in FIG. 1 and the asterisk denotes 

the complex conjugate of the quantity E (2) (r,m). Where 

E (2) ( r ,m) is the corresponding complex amplitude of the 
transitted beam. In the case of transparent media 50 
(Im(€^(Q)))=0): 



(6) 



55 



(7) 



for 5^ > 0 

/(r,m) tr - \ e • e 

far 5^ 0 

The dependence of transmitted light on the distance from the 
surface has been drastically changed as 5m cr becomes 

negative. At Sm^,^, the intensity I ( r ,n cr ) remains con- 
stant in infinite volume but begins to diminish with further 
decline of 8m cr . This volume becomes finite if the imaginary 65 
part of the dielectric constant is nonzero (Ime^co)^)), or, 
in other words, there is an absorbing species in solution. A 



finite (Ime^co) will to some extent smooth the singular 
behavior near 6m c pO. Nevertheless, if 



7m6 (2) (o))<li?e(e< 1 >(o»-€( l >(a)))l 



(8) 



60 



near the transmitted light thresholds, 8rn cr =0, the behavior of 
all reflected beams with different m will substantially 
change. Analytically, this behavior can be described by the 
following expressions for reflection coefficients R ( * ,p '(m), 
for s- and p-polarized light (transverse electric and trans- 
verse magnetic, respectively) in the interval of parameters in 
the threshold region (5m cr lies near zero): 



(9) 



20 



25 



Cf^(m) + 



CHm) N iS^I for d mcr * 0 
&{ p \m) \ I5 ffl I for 6,, >0 



as /me^-K) 

and for the phase reflected light: 



(10) 



S£ +(/meW) 2 



\ ( \l^ + tftHW + 8 ffl(y ) ) = 6?Xm) + di'Xm)^ 
as /me w -*0 

where parameters C li2t3A5 c,J,) (m) are constants and a iStP \n) 
corresponds to the phase of the reflected light For n*0, the 
constants C 1234 i5 (JJ,> (n) are related to the difference in 
dielectric constant between the grating lines (metal) and the 
substrate materials (glass). The important point is that the 
positions of singularities 8m cr =0 and the functional depen- 
dence on Sm cr in the interval 8m cr <l are dependent only on 
the bulk properties of the solution at distances from the 
solution-grating interface on the order of a few wavelengths. 
This means that this position, as well as the reflection 
coefficient's behavior in the vicinity of this position, in 
contrast to other properties of reflected light, will remain 
stable and unaffected by changing surface properties at 
distances less than 



II. Experimental 

Ttaiing to FIG. 1A, a single diffraction grating 101 was 
fabricated in a high-reflective chrome mask layer (1000 
angstrom thickness) deposited on a microscope slide 103. A 
Lloyd's mirror optical configuration was used to record 
interference fringes from an argon ion laser in a positive 
photoresist spin-coated on the slide. The resist was devel- 
oped to the chrome layer, and then the chrome was wet 
chemically etched to the glass substrate, forming a trans- 
mission grating. The grating period was measured to be 0.84 
microns with good uniformity and a sinusoidal surface 
profile. 
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Eleven ethanoi solutions 105 were prepared in deionized 
water from absolute ethanoi. The concentration range was 
\% EtOH by volume to 50% EtOH by volume. This 
corresponds to a refractive index range from 1 .33 1 to 1 .3579 
RI units. The refractive indices of the solutions were mea- 
sured with an Abbe refractometer 201 (Bausch and Lomb) to 
±0001 RI units at a wavelength of 589.3 nanometers and the 
dispersion was calculated, assuming normal dispersion in 
the visible region of the spectrum. 

The white light source 203 for these experiments was a 
tungsten bulb with a DC voltage source that was coupled 
into a 140 micron silica fiber optic 205. The fiber output was 
collimated via a 40X microscope objective and polarized 
with a Glan Taylor polarizing prism 207. The collimated 
light was incident on the backside of the substrate, and thus 
passed through the glass to the grating. Here the grating 
acted as the dispersing element and the m=-l reflected 
diffraction order was imaged on a 512 element fiber-optic 
window photodiode array 209 (EG&G Reticon). Approxi- 
mately 175 nm were imaged on the diode array, and the 
image was calibrated via. bandpass filters and exhibited a 
linear dependence of wavelength on diode number. The spot 
size of the beam was approximately 5 mm, and the diffracted 
light was focused on the detector with a 63 mm focal length 
achromat lens. The substrate itself was placed into a sand- 
wich flow cell which was mounted on a positioning stage 
that allowed translational adjustment and angular control to 25 
±0.083 degrees. 

As stated above, the substrate was mounted into a sand- 
wich flow cell and sealed with a silicone rubber gasket. The 
grating was in contact with the solution and the backside of 
the slide faced the optical train. The light that was collected 
on the photodiode array was reflected light, and thus did not 
pass through the solution. The optics were aligned such that 
a 37 degree coupling angle from normal was achieved. The 
polarization of the incident beam was chosen as transverse 
electric (TE), or transverse magnetic (TM). 

Measurements were taken in a stop flow format in order 
to remove flow effects. Sample was introduced into the cell 
with a peristaltic pump, and flushing sequence was a com- 
bination of well-defined flow and stop flow steps in order to 
allow complete flushing to minimize contamination from 
previous runs. Data was collected as full array images at 5 
second intervals, with 30 scans of the array taken for each 
solution. The solutions, numbered 1-1 1 , were run in random 
order with a water reference run between each alcohol 
solution. 

Each 30 scan block was placed into a Matlab™ (Math- 45 
works, Inc.) 211 file and averaged to obtain a single image 
for each solution run for calibration purposes. Each alcohol 
run was subtracted from its corresponding water blank, and 
this difference spectrum was the analytical signal that was 
used for Partial Least Squares (PLS) calibration. The indi- 
vidual spectra were then placed into a response matrix, and 
the refractive indices that were calculated from the Abbe 
refractomer dispersion measurements were used as the con- 
centration matrix for the PLS calibration. 

HI. Experimental Results 

FIGS. 2 and 3 are the experimental reflection spectra for 
m=l for TE and TM polarizations for the 11 solutions 
referenced to water. The plots show the dependence of the 
singularity position on dielectric constant of the solution. 
These spectra were loaded into the PLS algorithm, with 
measured refractive indices and dispersions used to calculate 60 
the refractive index for each solution at an average wave- 
length close to the singularity. For the TE polarization, the 
refractive index at 625 nm was used. For the TM polariza- 
tion, the refractive index at 614 nanometers was used for the 
calibration. The difference in the position of the singularity 65 
for the different polarizations arises from the adjustment of 
the optical train between the two polarization experiments. 



30 



35 
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Adjusting the polarizer effectively adjusts the incident angle. 
The calibration was also run using the alcohol concentra- 
tions for each solution in %vol/vol. EtOH in water. The 
calibration results are summarized in tables 1 and 2. The 
dependence of refractive index on EtOH concentration was 
calculated to a first approximation from a linear regression 
of the %EtOH and refractive indices at appropriate wave- 
lengths in order to get a final error of prediction in terms of 
%EtOH. It was found that the calibration was an order of 
magnitude better for the TM polarization and comparable 
for the TE polarization when the refractive indices were used 
in the calibration. This is probably due to the information 
contained in the feature in the TM plot at around 545 nm, 
which might be another diffraction order singularity or an as 
yet unexplained optical effect. The relationship between the 
refractive index and concentration of each solution yields a 
limiting error of 0.188% EtOH RMS due to the uncertainty 
in the Abbe refractometer of ±0001 RI units. The TE 
RMSECV shows that this limiting value was reached, so this 
sensor has achieved an accuracy that is as good as the 
standard method used to make up the calibration. The 
prediction plots are for the lowest factor with RMSEVC of 
less than 0.0001 RI units, as a lower error is not statically 
significant. The prediction plots are shown in FIGS. 4 and 5. 

TABLE 1 

Root mean square error of cross validation values for 
calibration of the GLRS sensor for EtOH solutions for 
TM incident polarization. 





RMSECV, 


RMSECV, 






% EtOH, 


RI units, 


Calculated 




calculated using 


calculated using 


RMSECV from RI 


# of PLS 


% EtOH as 


RIas 


calibration set, in 


Factors 


calibration set 


calibration set 


% EtOH units 


1 


1.3050 


0.000200 


0.379 


2 


1.0600 


0.000222 


0.421 


3 


0.6317 


0.000114 


0.216 


4 


0.2811 


0.000054 


0.102 


5 


0.3526 


0.000043 


0.082 


TABLE 2 


Root mean square error of cross validation values for 


calibration of the GLRS sensor for EtOH solutions for 




TE incident polarization. 






RMSECV, 


RMSECV, 






% EtOH, 


RI units, 


Calculated 




calculated using 


calculated using 


RMSECV from RI 


# of PLS 


% EtOH as 


RIas 


calibration set, in 


Factors 


calibration set 


calibration set 


% EtOH units 


1 


1.2605 


0.000585 


1.110 


2 


0.7944 


0.000343 


0.651 


3 


0.8105 


0.000376 


0.714 


4 


0.2427 


0.000106 


0.201 


5 


0.2286 


0.000097 


0.184 



The wavelength dependence of v To7I was determined 
using the dispersion data obtained from the Abbe refracto- 
meter and plotted for the 30% EtOH solution, TM polariza- 
tion, and plotted in FIG. 6. The important point to be made 
here is that the position of the singularity matches very 
closely with the minimum in the reflection coefficient plot, 
in the limit of experimental error and uncertainty. The 
detailed form of the reflection coefficient includes the con- 
stants mentioned in the introduction, but here we will not 
determine them as they must be deterrnined empirically for 
each solution. Thus, only the position of the minimum in the 
theoretical plot is of value for comparison. The theory 
developed here is with regard to the region close to the 
singularity, and any fitting to the actual experimental curve 
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at wavelengths removed from the vicinity of the singularity 
was not performed. The width of the minimum, according to 
Equation 9, is related to the imaginary part of the dielectric 
function. The theoretical plot assumes a value of zero for the 
imaginary part of the dielectric function. Thus, the shape of 
the theoretical curve at the singularity is much sharper than 
the experimental curve, as the experimental results indicate 
a non-zero dielectric function around the singularity. The 
offset of 7 nm between minima is due to experimental error 
in determining the grating period, coupling angle, and 
refractive index. This assumes that the grating lines are 
perpendicular to the plane of incidence. Deviation from this 
condition will result in spectral offset. Also, a one degree 
uncertainty in the coupling angle produces a 12 nm shift in 
the singularity position, and the difficulty in determining 
normal for the coupling angle reference is the most likely 
cause of the offset. However, the agreement in general form 
and position provides evidence that the theoretical predic- 
tions regarding the position of the singularity and the general 
form of the reflection spectrum are correct 

IV. Conclusion 

Even the relatively rough experiments demonstrated the 
feasibility of this proposed method. The calibration of the 
device via chemometrics indicates the usefulness of this 
method as a first order chemical sensor, yielding a vector of 
data for each sample. For analytical purposes it is possible 
to use the dependence of the singularity positions on con- 
centration as well as the form of the reflection and phase 
behavior near the singularities. If an absorbing species is 
present in solution in small amounts, the form of the 
reflection coefficient around the singularity will change. 
According to Equation 9 in the case of lme (2) =0 there must 
be a plateau in the dependence of RJ$ ma J on 8 m « near 
8 mBe =0. Thus the concentration of absorbing species can be 
measured according to the width of the pl ateau in R m (8 m J 

dependency near 8 moe =0 described by ^Ime (2) . By such 35 
methods it is possible, in principle, to measure relatively 
small concentrations of absorbing species. The obtained 
results show that the concentration sensitivity could be not 
less than Ime <2) =10" 2 . 

As this is not a surface plasmon technique, there is no 
materials constraint regarding the metal used for the grat- 
ings. Thus, gratings may be fabricated in any metal. This 
gives wide latitude for performing analyses, as the proper- 
ties of the metal may be chosen to impart selectively and 
exploit the dependence of the signal on dielectric differences 
between the metal and the substrate, thus making the grating 45 
an active chemical sensing region with which to perform the 
analyses. 

Also, the relative insensitivity to surface absorbents of the 
position of the singularity will allow bulk analyses to be 
performed in a medium that would not be amenable to an 
ATR orSPS probe. This bulk determination coupled with the 
reflection basis of the system makes it ideal for chromato- 
graphic applications where the pathlength is insufficient to 
perform transmission spectroscopic measurements, as is the 
case with microbore chromatography. 

As this is a reflection method, the analysis of average bulk 
properties in a highly scattering medium is possible. The 
signal would be dependent on the particle size and also on 
the homogeneity of the solution. This has application in 
processes or chemical systems that would not be suitable for 
transmission measurements. 

The relative ease of grating fabrication and also the ease 
with which this method can be coupled to a fiber optic 
allows for the integration of this sensor into a bulk optic that 
can be interrogated with an input fiber and an output fiber. 
Also, the relative signal is high due to a single pass reflection 
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measurement that does not undergo absorbance attenuation. 
Thus remote sensing with fiber optic coupling is achievable 
with this system. 

The importance advantage of the proposed method is its 
versatility. The method may be employed in the light fre- 
quency range spanning from the LTV to the mid IR, with 
wavelength selection determined by the grating period and 
angle of incidence. The potential here is for a single reflec- 
tion method with sensitivity in the UV and visible to allow 
for spectral measurements of absorbance and refractive 
index in liquids and gases in a variety of environments. 

While the preferred embodiment of the invention has been 
illustrated and described, it will be appreciated that various 
changes can be made therein without departing from the 
spirit and scope of the invention. 

The embodiments of the invention in which an exclusive 
property or privilege is claimed are defined as follows: 

1. An -analytical sensor for determining the index of 
refraction of a liquid or gas, said sensor comprising: 

a diffraction grating having a grating period formed on a 
dielectric and disposed in contact with said liquid or 
gas; 

a light source providing light incident on said diffraction 
grating at an angle of incidence and wherein such 
incident light does not traverse said liquid or gas; 

a photodiode array for capturing a reflection spectrum 
from said diffraction grating when said light is incident 
on said diffraction grating; and 

computing means for determining said refractive index 
based upon at least said angle of incidence, said grating 
period, and said reflection spectrum. 

2. The sensor of claim 1 wherein said light is polarized 
light that is transverse electric. 

3. The sensor of claim 1 wherein said light is polarized 
light that is transverse magnetic. 

4. The sensor of claim 1 wherein said index of refraction 
is a complex quantity having real and imaginary compo- 
nents. 

5. The sensor of claim 1 wherein said photodiode array 
can capture at least two diffracted orders of reflected light. 

6. The sensor of claim 1 wherein said liquid or gas 
includes solid particulates suspended therein as a matrix. 

7. A method for determining the index of refraction of a 
liquid or gas, said method comprising the steps of: 

placing in contact with said liquid or gas a diffraction 
grating having a grating period formed on a dielectric; 

illuminating said diffraction grating through the dielectric 
with light at an angle of incidence and wherein such 
incident light docs not traverse said liquid or gas; 

capturing a reflection spectrum from said diffraction grat- 
ing when said light is incident on said diffraction 
grating; and 

determining said refractive index based upon at least said 
angle of incidence, said grating period, and said reflec- 
tion spectrum. 

8. The method of claim 7 wherein said light is polarized 
light that is transverse electric. 

9. The method of claim 7 wherein said light is polarized 
light that is transverse magnetic. 

10. The method of claim 7 wherein said step of capturing 
said reflection spectrum captures at least two diffracted 
orders of reflected light 

11. The method of claim 7 wherein said liquid or gas 
includes solid particulates suspended therein as a matrix. 
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